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ABSTRACT 


A  model  of  double  diffusively  driven  interleaving  layers 
across  a  thermohaline,  barotropic  front  is  developed.  The 
mechanism  of  the  interleaving  Mas  proposed  by  Stern  (1967)  for 
small  perturbations  of  continuous,  large  scale  T-S  gradients  in 
the  saltfingering  regime  by  double  diffusive  interleaving.  The 
model  shows  that  this  mechanism  is  also  effective  for 
interleaving  across  an  initially  sharp  front,  after  strong  T-S 
inversions  devc^lop  between  intrusions,  in  regions  in  which  the 
mean  vertical  gradients  are  in  the  saltfingering,  diffusive,  or 
doubly  stable  regimes. 

The  vertical  scale  of  the  intrusions  increases  with 
increasing  cross-frontal  T-S  discontinuities,  and  decreases  with 
increasing  mean  vertical  density  gradients.  The  cross-frontal 
fluxes  caused  by  the  finescale  interleaving  can  be  substantial, 
and  the  vertical  finescale  fluxes  can  be  up  gradient. 
Sensitivity  studies  are  performed  to  study  the  effects  of  the 
mean  gradients,  microscale  flux  parameters,  and  intrusion 
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I. INTRODUCTION 

In  regions  with  large  scale  epipycnal  T-S  gradients  and  with 
diapycnal  gradients  suitable  for  salt  fingering,  the  salt 
fingering  drives  double  diffusive  finescale  interleaving  (DDI) 
which  undergoes  unstable  growth.  Stern  (1967)  developed  an 
analytic  perturbation  solution  for  this  phenomenon.  Advection 
along  the  large  scale  epipycnal  gradients  causes  alternate 
layers  of  ODI  to  become  relatively  warm  and  salty.  The 
thermohaline  properties  of  these  layers  cause  a  divergence  of  the 
vertical  flux  of  buoyancy  by  salt  fingering,  and  they  gain 
buoyancy.  If  the  layers  are  slightly  tilted  with  respect  to 
hori2ontal,  their  buoyancy  anomaly  will  lead  to  an  acceleration 
of  the  advection,  thus  closing  the  positive  feedback  loop  of 
unstable  growth. 

Many  observations  of  DDI  have  been  made  in  a  variety  of 
regions,  and  Stern's  (1967)  model  has  frequently  been  invoked  as 
a  possible  explanation  of  the  observations.  For  example,  Joyce 
et  al.  (1978)  observed  DDI  across  the  Antarctic  Polar  Front; 
Hal lock  (1985)  found  DDI  in  the  Norwegian  Sea;  and  Posmentier 
and  Houghton  (1978)  observed  DDI  across  the  shelf /slope  front  of 
the  Middle  Atlantic  Bight.  However,  there  were  two  significant 
gaps  between  Stern's  (1967)  explanation  and  these  observations: 
(1)  Stern's  explanation  applies  to  DDI  in  regions  of  continuous 
gradients,  not  necessarily  to  cross-frontal  DDI,  and  O)  Stern's 
explanation  applies  directly  only  to  the  initial  instability 
before  the  initial  vertical  gradients  become  significantly 
perturbed,  not  necessarily  to  the  later  stages  after  temperature 
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and  salinity  inversions  are  present. 

Toole  and  Georgi  (1981)  added  viscosity  to  Stern's  (1967) 
model,  and  found  the  solutions  to  be  qualitatively  unchanged 
although  the  DDI  developed  lees  rapidly.  They  also  noted,  as  did 
Posmentier  and  Hibbard  (1982),  that  an  "optimum"  tilt  of  the 
intrusive  layers  would  cause  a  maximum  growth  rate  of  the  DDI. 
Posmentier  and  Hibbard  (1982)  investigated  these  optimum  tilts 
and  found  that  the  tilt  component  parallel  to  the  mean  epipycnal 
gradients  had  relatively  little  effect  on  the  growth  rate  of  DDI, 
and  that  positive  growth  was  possible  with  along-gradients  tilts 

'  of  either  sign.  However,  the  cross-gradient  tilt  was  found  to  be 

critical  to  growth  rate. 

Posmentier  and  Hibbard  (1982)  also  found  that  very  large 
epipycnal  and  diapycnal  fluxes  of  salt,  heat,  and  buoyancy  are 
caused  by  the  DDI,  and  suggested  that  these  large  fluxes  might 
play  a  role  in  larger  scale  structures.  Posmentier  and  Kirwan 
(1985)  hypothesized  that  these  large  fluxes  caused  by  DDI  can  act 
to  enhance  the  temperature,  salinity,  and  buoyancy  signatures  of 
some  mesoscale  structures.  They  presented  quantitative  examples 
for  a  warm  and  a  cold  core  mesoscale  ring  to  demonstrate  the 
feasibility  of  the  hypothesis.  The  gap  between  Stern's  (1967) 
mechanism  and  the  frontal  regions  to  which  it  was  being  applied, 
however,  remained  an  unresolved  issue. 

McDougall  (198Sa,  198Sb)  suggested  a  model  which  moves 

toward  removing  the  limitations  of  the  perturbation  assumptions. 

I 

The  model  he  proposed  is  an  alternating  slabs  model,  rather  than 
the  continuous  harmonic  function  of  space  upon  which  the  earlier 
work  had  been  based.  It  was  demonstrated  that  if  the  temperature 
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and  salinity  discontinuities  between  successive  slabs  are 
initially  all  equal,  the  system  will  be  unstable  with  respect  to 
temperature  and  salinity  perturbations  of  opposite  signs  in 
alternate  layers.  The  growth  mechanism  is  identical  to  that 
originally  proposed  by  Stern  (1967).  However,  these 
perturbations  can  continue  to  grow  beyond  small  amplitudes  until 
temperature  or  salinity  inversions  occur.  Furthermore,  it  was 
shown  that  a  steady  state  solution  exists  for  DDI  after  the 
inversions  have  become  established.  However,  this  model  was  no 
more  applicable  to  fronts  than  was  Stern's  (1967)  earlier  model. 

Turner  (1978)  proposed  a  mechanism  similar  to  Stern's  (1967) 
to  explain  both  laboratory  and  oceanic  observations  of  cross- 
frontal  DDI.  This  was  elaborated  upon  by  Ruddick  and  Turner 
(1979),  who  also  explained  the  dependence  of  the  vertical  scale 
of  the  DDI  on  the  cross-frontal  thermohaline  gradients  and  on  the 
vertical  density  gradient.  In  these  models,  however,  the  earth's 
rotation  was  neglected,  so  their  results  can  be  applied  to 
oceanographic  observations  only  by  extrapolation. 

In  this  paper,  a  model  of  double  diffusively  driven 
interleaving  across  a  vertical,  thermohaline,  barotropic  front  is 
formulated.  One  purpose  of  the  model  is  to  demonstrate  that  the 
original  mechanism  proposed  by  Stern  (1967)  can  produce  such 
DDI,  even  though  the  original  theory  and  its  extensions  were  not 
directly  applicable  either  to  fronts  or  to  interleaving  in  which 
the  vertical  salinity  and  temperature  gradients  were  perturbed 
beyond  small  amplitudes.  The  second  purpose  is  to  demonstrate 
that  all  these  solutions  exist  not  only  for  large  scale 
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environments  which  are  conducive  to  salt  fingering,  but  also  in 
regions  with  T  and  S  both  decreasing  upward,  and  in  doubly  stable 
regions.  Third,  the  model  will  be  used  to  evaluate  several 
observable  and  significant  properties  of  cross-frontal  DDI ,  such 
as  horizontal  length  scales,  strength  of  the  salinity  and 
temperature  variations,  and  cross-frontal  fluxes  caused  by  the 
DDI.  Finally,  the  model  will  be  used  to  evaluate  the  sensitivity 
of  cross-frontal  DDI  to  the  hydrographic  features  of  the  front, 
to  the  parameters  controlling  vertical  microscale  (saltfingering 
and  diffusive)  fluxes,  and  to  the  vertical  thickness  and  tilts  of 
the  layers. 

II.  FORMULATION  OF  THE  MODEL 
A.  Structure  and  Form 

The  front  and  its  environment,  before  the  front  is  disturbed 
by  intrusions,  consists  of  two  vertically-stratified  regions  with 
the  same  vertical  gradients  of  salinity,  temperature  and  density, 
Sz»  ^z»  ^*z  *  respectively.  The  vertical  front  is  in  the  plane 

x=0  in  the  x(east),  y<north),  z<down)  coordinate  system.  There 

is  a  discontinuity  in  both  salinity  and  temperature,  of  5^  and 
Tj,  respectively,  crossing  the  front  from  the  x<0  side  to  the  x>0 
side,  but  no  discontinuity  in  density.  The  ambient  conditions 
are  described  by  S^Sq-Si+S^z  ^  TsTq-Tj+T^z,  and  p=Po'*’Pz2 
x<0,  the  cold  fresh  <CF)  side  of  the  front,  and  by  S=Sq-Si+S2Z, 
T^Tq-Tj+T^z,  and  p=Po'*'9z*  x>0,  the  warm  salty  <WS>  side  of 

the  front. 

The  front's  shape  is  changed  by  a  series  of  equally  thick, 
parallel,  tilted  slabs  of  water  moving  in  opposite  directions 
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(Fig.  1).  The  even-numbered  slabs  carry  CF  water  a  distance  SL 
in  the  positive  x-direction,  and  the  odd-numbered  slabs  bring  US 
water  in  the  opposite  direction.  Microscale  double  diffusive 
fluxes  of  salt  and  heat  between  adjacent  layers  release  potential 
energy  and  change  the  densities  of  the  layers,  causing  pressure 
gradients  which  drive  the  finescale  interleaving  motion.  The 
layers  are  separated  by  planes  of  constant  kx+hy+z,  and  the  value 
of  the  plane's  constant  is  (j-'4)H  at  the  top  of  the  j  '  th  layer, 
where  H  represents  the  layers'  thickness. 

The  X-,  y-,  and  z-components  of  velocity  in  the  CF  layers 
are  U,  V,  and  W.  In  the  WS  layers,  they  are  -U,  -V,  and  -W.  One 
condition  for  the  validity  of  the  solution  for  velocity  is  that  U 
must  be  positive.  The  salinity,  temperature,  and  density  in  the 
CF  n'th  layer  (in  this  report,  n  is  always  used  to  denote  an  even 


integer)  are 

S=Sq— Sj  ■►nHSj'^Sj^  x+Sy y— S  '  (1) 

T=TQ-Ti+nHT2+T^x+Tyy-T'  (S) 

p=P0'*'"^^<?z'^9x’<'*'Pyy“«''  ^3) 

In  the  n+l'st  layer  (which  is  WS) ,  they  are 

S— Sq+S  j +nHS2'*"S,<x+Syy+S '  (4) 

TsTo+Ti+nHTj+T^x+Tyy+T'  (5) 


p=PO'*’"'^Pz‘^PxX'^Pyy*'P’  (6) 

A  linear  equation  of  state  is  used,  so 

Px=(iS^-0(T^  (7) 

Py=flSy-«Ty  (8) 

P2=<1S2-<xT2  (9) 

p'=(js'-«T'  no) 


6 


0=(lS  j  ~<xT  j 


(11) 


where  (I  is  the  coefficient  of  haline  contraction,  and  «  is  the 
coefficient  of  thermal  expansion. 

B.  Boundary  Conditions 

The  water  entering  a  CF  layer  at  x=-L  is  water  from  the 
adjacent  ambient  environment  in  the  vertical  range  of  2H  centered 
a  distance  H'  below  the  center  of  the  CF  layer.  Similarly, 
the  water  entering  a  MS  layer  at  x=-<-L  is  water  from  the  adjacent 
ambient  environment  in  the  vertical  range  of  2H  centered  a 
distance  H'  above  the  center  of  the  MS  layer.  The  process  at  the 
origin  of  the  CF  layers,  at  x=-L,  leads  to  the  following 
boundary  conditions: 

-S„L-S'=S2kL+S2H' ,  and  Sy=-hS2  (12,13) 

-T^L-T'sTjkL+TjH* ,  and  Ty=-hT2  (14,15) 

-p,<L-p  '  9y“”^^P2  (16,17) 

Eqs  (16,17),  however,  are  not  independent.  They  follow  from  the 
four  prior  boundary  condition  equations  and  the  equation  of 

state.  The  boundary  conditions  arising  from  the  process  of  MS 
layer  formation  at  x^-i-L  are  redundant  with  those  just  given. 


The  differences  in  S,  T,  and  p  across 

the  interfaces 

separating  each  CF  layer 

from  the  MS  layer  above 

it,  or  below  it, 

respectively  are 

Sp_i -Sn=2Si -HSj+ES ' 

and  S^+ j-S,^=2Si+HS2+2S  ' 

( 18, 19) 

Tn-l-Tn=eTi-HT2+2T- 

and  T„+i-Tn=2Ti+HT2+2T' 

(20,21 ) 

Pn- 1 ~  Pn“2pi -Hp2+2p • 

and  pp+i-p„=2pi+Hp2+ap‘ 

(22,23) 

Another  condition  for  the  validity  of  the  results  is  that  all  the 
discontinuities  in  Eqs  (18)-(21)  must  be  positive,  because  the 
microscale  flux  laws  to  be  introduced  below  are  predicated  on 
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this  condition 


Furthermore,  stability  requires  that  the  density 


discontinuity  in  Eq  (SS)  be  negative,  and  in  Eq  (S3)  positive. 

It  is  necessary  that  the  vertical  average  of  density  over  a 
range  SH  should  be  the  same  immediately  to  the  left  of,  and 
immediately  to  the  right  of  x=+L.  If  this  were  not  true,  a 
horizontal  pressure  discontinuity  would  exist  at  x=+L  which  would 
increase  ad  inf ini  turn  with  depth.  At  (x,y)  =  (L,y),  the  average 
density  inside  the  n  ‘  th  and  n+l‘st  layer  is  pq+ ( n+'4 )  Hp2'*'PxL+Pyy . 
In  the  ambient  water  on  the  other  side  of  the  discontinuities  at 
the  layer  ends  and  x^l,  the  average  density  over  the  same  depth 
range  is  •  Equating  these  two  average 

densities  results  in  Eq  (17)  again,  and  in  the  new  result 

p^=-kp2  (24) 

A  similar  requirement  at  x=-L  results  in  the  two  identical 
equations. 

While  the  vertically-averaged  pressure  discontinuity  across 

the  planes  x^-i-L  and  x=-L  must  be  zero,  the  discontinuity  across 

these  planes  is  not  zero  if  averaged  only  over  those  layers  which 

are  entraining  ambient  water  and  moving  it  toward  x=0.  The 

average  pressure  discontinuity  at  xs-l  at  the  depths  of  the  CF 

layers  (and  at  x^-^-L  at  the  WS  layers),  according  to  Bernoulli's 

2  2 

Law,  must  be  p(U  +V  ) .  It  can  be  shown  from  the  hydrostatic  law 
that  the  pressure  in  the  ambient  fluid  is 

P=PQ+PQg  (  z-zo)‘*’P2g  <  (25) 

where  zq  is  the  depth  of  the  center  of  the  n'th  (CF)  layer  at 
x=-L,  and  Pq  is  the  pressure  at  that  point.  Now,  if  the  ambient 
water  in  the  depth  range  Zq+H'-H  to  zq+H  '  •♦•H  is  vertically  mixed 
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before  entering  the  CF  layer,  the  pressure  in  this  vertically- 
mixed  ambient  fluid,  in  which  the  density  is  Po'*'9z  ^  ^  »  will 
be 


P=Po+Po9<2"='o'‘^P29‘^2<zo+H'  )  (z-zo)-(H+H'  )  (H-H‘  )  3/2  (26) 


if  the 

pressure 

at 

the 

too  (  i .e. . 

at  zq+H'-H)  is  the 

same  as 

before 

mixing. 

On 

the 

other  hand. 

the  pressure  inside 

the  CF 

layer,  which  has  the  same  density,  will  be 

P=Po+Pzg^2<Zo+H' ) (z-Zo)+<H+H' /4)H) 3/2  (27) 

if  the  pressure  at  the  top  is  the  same  as  in  the  ambient  water. 
The  difference  between  the  pressures  in  Eqs  (26)  and  (27), 
averaged  over  the  depth  range  of  the  CF  layer  from  zq-H/2  to 
zo'‘‘H/2,  is 

P=P29  <  3H+2H • ) ( H-2H ■  ) /S  ( 28 ) 

According  to  Bernoulli's  law,  this  pressure  difference  can  be 
equated  with  the  increase  in  kinetic  energy  per  unit  volume  when 
the  ambient  water  at  rest  acquires  the  velocity  of  the  CF  layer 
upon  entering  the  layer.  Thus, 

P2g(3H+2H'  )(H-2H'  ) /8=C  (  Pq+Pz  <  '  )3(U^+V^)  (29) 

Neglecting  the  smaller  part  of  the  density  on  the  right  side  of 
the  Eq  (29)  results  in  the  quadratic  equation  for  H', 

H ■ '  +  ( Po  <  ) / p2g-3H^/4 ) =0  ( 30 ) 

The  same  equation  follows  from  the  application  of  Bernoulli's  Law 
to  the  pressure  discontinuity  at  x=-t-L  at  the  depths  of  the  MS 
1 ayers . 

H'  may  also  be  related  directly  to  the  density  anomaly  in 
the  US  layers  by  combining  Eqs  (16)  and  (24)  to  obtain 

p's-PjH'  (31) 
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C.  Equations  of  Motion 

In  anticipation  of  applying  Newton's  Second  Law  to  the 
dynamics  of  the  layers,  the  pressure  gradients  within  the  layers 
can  be  derived  from  the  hydrostatic  law  and  the  assumption  that 
the  gradients  are  equal  and  opposite  in  the  n'th  (even,  CF )  and 
n+l‘st  (odd,  WS)  layers: 

T8p78xTp=-gkp‘  and  T8P7dyTp=-ghp *  (33,33) 

( 8P7dx7p+i=gkp ■  and  < dP7dyTp+i=ghp '  (34,35) 

where  g  is  the  acceleration  of  gravity.  The  equations  of  motion 
in  the  x-  and  y-directions  respectively  for  the  n'th  layer,  which 
are  redundant  with  those  for  the  n+l'st  layer,  are 

gkp' +(-A/H)U+pfV=0  (36) 

ghp'  +  (-*A/H)V+<-pf  )U=0  (37) 

assuming  that  U  and  V  are  constants,  where  A  is  the  viscosity  and 
f  is  the  Coriolis  parameter.  If  the  velocities  are  parallel  to 
the  planes  separating  the  layers,  the  z-component  of  velocity 
must  satisfy 

W=-(kU-»-hV)  (38) 

D .  Conservation  Equations 

The  microscale  flux  laws  governing  the  vertical  fluxes  of 
salt  and  heat  between  layers  are  required  for  the  conservation 
equations.  Assuming  that  saltfingering  is  the  cause  of 

microscale  fluxes  from  the  n-1 ' th  (odd,  WS)  layer  down  to  the 
n'th  layer,  the  downward  salt  and  heat  fluxes  across  the 
interface  is 

Fg=pK(Sp_i-Sp)  and  Fj=t  ((5/a)  pK(Sp_i-Sp)  (39,40) 

where  K  is  the  coefficient  of  diffusion  for  the  salt  fingering 
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interfaces,  and  t  is  the  ratio  of  density-equivalent  heat  to  salt 
carried  down  by  the  salt  fingering.  If  the  interface  below  the 
n'th  layer  is  a  diffusive  interface,  the  upward  fluxes  across 
this  interface  from  the  n+l'st  layer  below  will  be 

Fg=T ' <a/a) pK* (Tp+i-Tp)  and  Fj=pK ' ( Tp+i-Tp )  (41,48) 

If  S  and  T  are  steady-state,  the  equations  for  conservation  of 
salt  and  heat,  respectively,  in  the  n’th  layer,  after 

substituting  Eqs  (IS-Sl)  into  (39-42)  and  adding  advection  terms, 
are 

(K/H>  (2Si+2S’-HSz  )  +  («/(!  )(t’K’  /H)  (  2Ti+2T  ’ +HT2  ) 

-US,<-VSy=0  (43) 

(fl/a) (tK/H) (aSi+2S’-HS2)+<K’/H) ( aT^ +2T ’ +HT2 > 

-UTx-VTy=0  (44) 

which  are  redundant  with  the  equations  for  the  n+l'st  layer. 

III.  SOLUTION  OF  GOVERNING  EQUATIONS 

In  the  preceding  section,  twelve  (12)  constants  and  more  or 
less  arbitrary  parameters  have  been  used: 

H  Si  Sj  Ti  T2  K  T  K'  T'  A  k  h 
After  substituting  Eqs  (12-15)  into  (43,44),  there  remain  eight 
(B)  unknowns: 

S'  T'  p'  U  V  W  L  H* 

which  appear  in  eight  (8)  independent,  nonlinear  equations: 

10  30  31  36  37  38  43  44 

The  following  pseudocode  represents  a  procedure  which  may  be  used 
to  solve  the  equations  numerically. 

Assign  values  to  constants  and  parameters 
Guess  values  of  U  and  V 


Solve  (30)  for  H 


Do  until  H'  stops  changing; 

Find  p ‘  from  < 31 ) 

Find  U,  V,  U  from  (36-38) 

Solve  (30)  for  H‘;  let  new  H'=average  of  old  and  new 
End  do  loop 
Find  p '  from  ( 31 ) 

Guess  lower  and  upper  bounds  for  L 

Find  S‘  and  T‘  from  (43,44)  at  each  bound  of  L 

Find  corresponding  p‘  values  from  (10) 

Do  until  L  stops  changing: 

Interpolate  to  new  L  at  which  p'  agrees  with  (31) 

Find  p*  at  new  L  using  (43,44)  and  then  (10) 

Replace  one  of  bounds  for  L  by  new  L 
End  do  loop 

Convergence  occurs  in  only  a  few  steps.  However,  for  some  values 
of  the  constants  and  parameters,  convergence  does  not  occur 
within  reasonable  bounds  for  L  (for  example,  between  1  cm  and 
1000  km),  or  convergence  occurs,  but  one  or  more  of  the 
assumptions  mentioned  in  the  preceding  section  (for  example,  U>0) 
is  violated. 

IV.  NUMERICAL  EXPERIMENTS 

A  number  of  numerical  experiments  were  performed  to 
demonstrate  the  existence  of  solutions  for  realistic  constants 
and  parameters,  and  to  find  the  sensitivity  of  the  solutions  to 
the  constants  and  parameters.  The  assumed  values  together  with 
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results  describing  the  solutions  for  11  different  experiments  are 
summarized  in  Table  1. 

Table  1.  includes  among  the  results  the  values  of  the 

salinity  and  temperature  changes  across  the  saltfingering 

interfaces  between  each  US  layer  and  the  CF  layer  immediately 

below  it  < aS+  and  aT+,  respectively),  and  of  the  diffusive 

interfaces  between  each  WS  layer  and  the  CF  layer  above  < aS-  and 

aT-).  The  magnitudes  of  these  four  discontinuities  denote  the 

"strength"  of  the  DDI  from  the  perspective  of  hydrographic  data 

interpretation.  Also  included  are  the  velocity  components  (U,  V, 

W),  the  half-length  of  the  intrusion  <L) ,  and  the  vertical 

displacement  of  ambient  water  being  drawn  into  the  interleaving 

layers  <H').  In  addition,  the  cross-frontal,  horizontal  salt 
2  2 

fluxes  (mg/cm  -s)  and  heat  fluxes  (cal/cm  -s)  caused  by  the  DDI, 

and  the  effective  horizontal  diffusion  coefficients  for  salt  and 

2 

heat,  K"S  and  K"T  (cm  /s),  are  included.  The  diffusion 
coefficients  are  defined  here  as  the  fluxes  divided  by  the 
horizontal  gradients  in  the  frontal  zone,  Sj/L  and  Tj/L, 

respectively. 

In  all  the  numerical  experiments,  the  constants: 

«  (1  So  Tq  g  f 

were  set  to 

1.4E-4  7E-4  35  15  980  lE-4 

respectively,  in  cgs  units  except  for  salinity,  which  is  in  mg/g. 

Experiment  A  (see  Table  1,  column  A)  is  used  as  a  reference 
for  comparison  with  the  other  experiments  reported  here.  In  this 
experiment,  a  front  with  a  Img/g  salinity  discontinuity  and  a  5®C 
temperature  discontinuity  separates  two  regions  with  temperature 
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and  salinity  decreasing  downward  —  saltfingering  regions.  The 
1000  cm  thickness  of  the  layers  is  typical  of  observed  values 
(Posmentier  and  Houghton,  1978).  The  tilts  are  consistent  with 
those  discussed  by  Posmentier  and  Hibbard  (198S>.  The  values  of 
the  process-controlling  parameters  are  representative  of  those 
generally  discussed  in  the  literature,  if  the  interfaces  between 
layers  are  100  cm  thick. 

The  results  indicate  that  a  solution  exists:  the  mechanism 
under  discussion  can  cause  DDI  under  the  conditions  of  experiment 
A,  which  are  conditions  not  unlike  those  observed  in  various 
frontal  regions.  Furthermore,  the  results  are  comparable  to 
DDI  observed  in  such  areas  -  the  large  S  and  T  discontinuities 
across  the  saltfingering  interfaces  below  the  US  layers,  and 
across  the  weaker  diffusive  interfaces  above  the  US  layers,  are 
detectable,  and  the  total  length  of  the  layers  (lS.6Skm)  makes 
their  observation  certain  in  any  frontal  field  study. 

The  cross-frontal  fluxes  calculated  by  the  model  are 
substantial,  and  indicate  that  cross-frontal  DDI  may  be  an 
important  factor  in  transport  between  water  masses  bounded  by 
fronts.  In  addition,  the  vertical  fluxes  and  salt  and  heat  are 
upward,  which  is  the  upgradient  direction.  The  effective 
coefficients  of  vertical  diffusion  are  negative*.  The  model  thus 
lends  qualitative  support  to  the  underlying  assumption  of 
Posmentier  and  Kirwan's  counterintuitive  hypothesis  (1985)  that 
DDI  can  amplify  the  contrast  between  contiguous  water  masses. 

The  advection  of  heat  and  salt  into  our  out  of  each  layer  is 
exactly  balanced  by  the  microscale  diffusion  out  of  or  into  the 
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layer,  as  required  by  Eqs  (43,44).  In  experiment  A  and  most  of 
the  other  experiments,  the  pressure 

gradient  acceleration  in  the  CF  layers  is  oriented  somewhat 
closer  to  the  negative  x-direction  (out  of  the  frontal  region) 
tt~  a-  the  y-direction  (along  front)  because  of  the  assumed 
components  of  tilt.  The  flow  is  roughly  in  geostrophic 
equilibrium,  but  viscosity  is  not  negligible,  so  the  velocity  is 
aligned  somewhat  closer  to  the  y-direction  (along  front)  than  the 
x-direction  (cross  front). 

An  additional  set  of  experiments  was  performed  in  which  the 
value  of  S2  was  varied,  keeping  82  and  T2  both  negative  and  the 
ratio  S2/T2  constant,  so  conditions  remained  favorable  for 
saltfingering.  In  these  experiments,  DDI  occurred  only  for  -S2>- 
2.9E-S  mg/g/cm.  As  the  magnitudes  of  the  vertical  gradients  were 
increased  beyond  this  minimum,  the  strength  of  the  interleaving 
increased  relatively  faster  than  the  gradients,  and  L  decreased 
less  rapidly  than  the  gradients'  increases.  Experiment  B  is  an 
example  of  these  effects.  (See  Table  1.)  A  further  consequence 
of  the  increased  vertical  gradients  in  experiment  B  is  the  change 
in  the  range  of  H  for  which  DDI  occurred:  100<H<1390  cm.  The 
lower  limit  of  H  decreased  considerably  compared  to  experiment  A, 
while  the  upper  limit  decreased  very  slightly. 

The  magnitudes  of  the  cross-frontal  discontinuities  affect 
the  results  in  the  same  sense  as  the  vertical  gradients,  but  not 
as  strongly.  (See  experiment  C  in  Table  1.)  Interleaving  occurs 
under  the  same  conditions  as  experiment  A  but  with  different  S^ 
and  Tj  (satisfying  Eq  (11),  which  states  that  the  front  is 
barotropic)  only  if  Si>0.30.  However,  doubling  S^  causes  only 
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small  increases  in  the  strength  of  the  stronger  saltfingering 
interfaces,  but  causes  near  doubling  of  the  salinity 
discontinuity  and  more  than  doubling  of  the  temperature 
discontinuity  in  the  Meaker  diffusive  interfaces  above  the  WS 
layers,  making  detection  more  likely. 

A  further  consequence  of  the  increased  cross-frontal 
discontinuities  in  experiment  C  is  that  L  nearly  doubles,  making 
it  easier  to  find  the  intrusions.  In  addition,  the  range  of  H 
for  which  DDI  occurs  changes  to  S70<H<2990  cm.  The  lower  limit 
of  H  decreased  very  slightly  compared  to  experiment  A,  while  the 
upper  limit  doubled. 

Changing  the  layer  thickness  from  experiment  A  precluded 
DDI  unless  2900<H<1460  cm.  Thicker  layers  (experiment  D)  tend  to 
cause  somewhat  stronger  saltfingering  interfaces  and  longer 
intrusions,  but  cause  the  weaker  diffusive  interfaces  to  approach 
disappearance.  It  is  particularly  interesting  that  there  is  a 
definite  range  of  values  of  H  outside  of  which  interleaving 
cannot  occur,  a  phenomenon  which  is  implicit  in  the  literature 
but  has  not  been  explained  previously. 

The  cross-frontal  tilt  of  the  layers,  k,  must  satisfy 
0.0094< Jk |<0.0250  for  DDI  to  be  possible  without  changing  other 
parameters  in  experiment  A.  Doubling  the  magnitude  of  k 
(experiments  E  and  F)  causes  a  ninefold  increase  in  iS-  and  an 
eighteenfold  increase  in  iT-,  which  are  the  least  discernible  of 
the  discontinuities,  thus  making  the  interleaving  much  more 
noticeable  in  hydrographic  sections  or  T-S  diagrams.  However, 
the  doubling  causes  2«L  to  decrease  to  3.8  km,  considerably 
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smaller  than  2*L  in  experiment  A,  but  still  not  difficult  to 
locate  in  the  field. 

Changing  the  sign  of  k  (experiment  F)  has  the  drastic  effect 
of  reducing  U  and  L  nearly  to  zero,  making  the  intrusions  too 
slow  to  be  of  any  significance  in  causing  cross-frontal  fluxes 
and  too  slow  to  find  in  the  field.  This  occurs  because,  if  k>0, 
the  cross-frontal  velocity  of  the  anomalously  low  density  water 
entering  the  US  layers  would  cause  them  to  sink.  In  order  to 
rise,  these  layers  must  move  predominantly  in  the  y-direction, 
and  are  close  to  geostrophic.  This  defeats  the  basic  mechanism 
which  causes  the  interleaving,  causing  the  layers  to  be  only  four 
times  longer  than  their  thickness. 

The  along-frontal  tilt,  h,  is  quite  critical  to  the 
interleaving.  If  the  value  of  h  in  experiment  A  is  varied, 
DDI  is  still  possible  for  O.OO10<-h<O.O23O  <h  must  remain 
negative.)  The  most  significant  effect  of  a  larger  tilt  (see 
experiment  G  in  Table  1)  is  the  disproportionate  increase  in  iS- 
and  in  iT-,  the  least  discernible  of  the  discontinuities,  thus 
having  the  same  effect  as  increasing  the  magnitude  of  k 
(experiment  E)  in  making  the  interleaving  much  more  noticeable  in 
hydrographic  sections  or  T-S  diagrams.  In  addition,  the  flow  is 
closer  to  geostrophic.  Unlike  increases  in  k,  however, 
increasing  h  causes  only  a  small  decrease  in  5*L. 

DDI  continues  to  occur  if  K  in  experiment  A  is  varied  in  the 

p 

range  0<K<.01  cm  /s.  Reducing  K  (see  experiment  H)  has  the 
effect  of  increasing  the  strength  of  the  intrusions,  especially 
of  the  weak  discontinuities  across  the  diffusive  interfaces  above 
the  US  layers.  L,  however,  is  reduced,  but  remains  large  enough 
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for  the  intrusions  to  be  found  in  the  field.  The  upper  limit  on 
K  in  the  series  of  experiments  including  A  and  H  should  not  be 
generalized,  because  all  other  parameters  including  k  and  h  were 
kept  constant  in  this  particular  sensitivity  study,  whereas  in 
reality  new  values  of  k  and  h  can  permit  DDl  with  K  larger  than 
the  upper  limit  in  this  series  of  experiments. 

In  experiment  1,  t  was  increased  from  experiment  A.  Other, 
similar  experiments  found  DDI  to  occur  for  0.51<t<1.0.  Larger 
values  of  t  have  the  same  effect  on  the  intrusions  as  smaller 
values  of  K.  The  action  common  to  both  increasing  t  and 
decreasing  K  is  that  of  decreasing  the  rate  at  which  density  in 
the  US  layers  is  lowered  by  saltfingering  downward  into  the  CF 
layer  below.  This  is  the  reason  for  their  common  effect  on 
increasing  the  intrusion  strength  while  decreasing  intrusion 
length . 

Changing  the  viscosity  in  experiment  A  within  the  range 
0.019<A<0.8  cm/s  permitted  DDI  to  occur.  Experiment  J  shows  that 
larger  A  causes  the  interleaving  velocity  to  decrease  only 
slightly  in  speed,  while  changing  in  direction  towards  the  x- 
axis.  This  is  because  H'  is  required  to  increase,  increasing  A' , 
which  in  turn  increases  the  pressure  gradients  within  the  layers, 
offsetting  most  of  the  tendency  to  slow  with  increasing 
viscosity.  At  the  same  time,  the  velocity  veers  further  from 
geostrophic  towards  counter  gradient.  Paradoxically,  the 
interleaving  strength  increases  (also  because  of  the  increase  in 
H'),  and  S*L  is  reduced. 

In  every  experiment  described  so  far,  the  ambient  water 
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outside  the  interleaving  region  was  in  the  saltfingering  regime. 
Although  this  is  a  necessary  condition  for  DDI  in  Stern’s  (1967) 
model  for  interleaving  within  a  region  with  continuous  gradients, 
it  is  not  a  necessary  condition  for  cross-frontal  interleaving. 
Ambient  vertical  gradients  in  the  diffusive  regime  or  in  the 
doubly  stable  regime  are  also  consistent  with  cross-frontal 
interleaving  (the  fourth  possibility,  doubly  unstable,  is  of  no 
interest.)  In  fact,  experiments  in  these  regimes  suggest  that 
interleaving  is  stronger  and  can  occur  over  a  wider  range  of  the 
other  parameters  (Sj^,  h,  K,  etc.),  than  in  the  saltfingering 
regime.  One  example,  experiment  K,  is  included  in  Table  1.  In 
this  experiment,  the  magnitude  (but  not  the  sign)  of  T2  and  the 
value  of  ^2  are  the  same  as  in  experiment  A,  but  82  is  changed. 
The  interleaving  is  stronger  than  in  any  experiment  except  H,  and 
S«L,  although  smaller  than  in  experiment  A,  is  still  not  an 
impediment  to  finding  intrusions  in  a  field  study. 

V.  CONCLUSIONS 

1.  A  model  extending  the  mechanism  proposed  by  Stern  (1967)  is 
capable  of  explaining  interleaving  across  thermohaline  fronts, 
with  alternating  US  and  CF  layers  separated  by  saltfingering  and 
diffusive  interfaces.  The  model's  results  and  the  observed 
features  of  cross-frontal  DDI  are  quantitatively  comparable. 

2.  Concl  ision  #1  applies  to  sal tf inger ing  regimes  (regions  with 
salinity  and  temperature  gradients  both  directed  upward),  to 
diffusive  regimes,  and  to  doubly  stable  regimes. 

3.  The  cross-frontal  fluxes  produced  by  DDI  can  be  significant 
to  salt  and  heat  budgets  of  the  water  masses  separated  by  the 
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front.  The  vertical  components  of  the  fluxes  are  up-gradient ( 

4.  The  upper  bound  of  the  range  of  possible  vertical  scales 
increases  approximately  in  proportion  to  the  cross-frontal 
salinity  and  temperature  discontinuities,  and  decreases  at  the 
lower  bound  somewhat  more  slowly  than  the  inverse  of  the  vertical 
density  gradient  This  is  qualitatively  consistent  with  the 
results  of  Ruddick  and  Turner  (1979). 

5.  The  strength  of  the  DDI,  as  manifested  by  the  magnitudes  of 
the  salinity  and  temperature  discontinuities  between  layers,  and 
the  horizontal  length  scale,  are  sensitive  to  the  vertical 
gradients,  to  the  cross-frontal  discontinuities,  to  the  layer 
thickness  and  both  components  of  tilt,  and  to  the  parameters 
controlling  the  saltfingering  and  double  diffusion  across  the 
interfaces  between  layers. 
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TABLE  1.  SUMMARY  OF  NUMERICAL  EXPERIMENTS 


A 

8 

C 

D 

E 

F 

H(cm) 

1000 

1000 

1000 

1400 

1000 

1000 

SI 

0.5 

0.5 

1 

0.5 

0.5 

0.5 

T1 ( xC) 

2.5 

2.5 

5 

2.5 

2.5 

2.5 

Sz ( 1/cm) 

-4E-05 

-1 .6E-04 

-4E-05 

-4E-05 

-4E-05 

-4E-05 

T  z ( xC/cm ) 

-0.00068 

-0.00272 

-0.00068 

-0.00068 

-0.00068 

-0.00068 

K ( cm/s ) 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

g 

0.56 

0.56 

0.56 

0.56 

0.56 

0.56 

K ' (cm/s) 

0.002 

0.002 

0.002 

0.002 

0.002 

0,002 

g' 

0.56 

0.56 

0.56 

0.56 

0.56 

0.56 

A(cm/s) 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

k 

-0.01 

-0.01 

-0.01 

-0.01 

-0.02 

0.02 

h 

-0.004 

-0.004 

-0.004 

-0.004 

-0.004 

-0.004 

aS+ 

0.103 

0.552 

0.  119 

0.132 

0.281 

0.281 

aT  + ( xC ) 

1  .408 

5.947 

1.490 

1  .91 1 

2.209 

2.210 

aS- 

0.023 

0.232 

0.039 

0.020 

0.201 

0.201 

aT-( xC) 

0.048 

0.507 

0.130 

0.007 

0.849 

0.849 

U(cm/s) 

1  .59 

4.87 

1  .59 

2.05 

1  .66 

0.018 

V(cm/s) 

-2.46 

-7.53 

-2.46 

-3.59 

-4.04 

4.37 

W(cm/s) 

0.0061 

0.0185 

0.0061 

0 . 0062 

0.0171 

0.0171 

L(  km) 

6.31 

2 . 62 

10.87 

8.91 

1  .9 

0.02 

H’ (cm) 

430.8 

329.3 

430.8 

601  .8 

338.9 

338.9 

S  Flux 

0.0501 

0.9545 

0.0628 

0 . 0779 

0.2000 

0.0022 

K"S 

6 . 32E+04 

5.00E+05 

6.83E+04 

1 .39E+05 

7.60E+04 

8.68E+00 

T  Flux 

0 . 5788 

7.8577 

0.6440 

0.9830 

1 .2691 

0.0138 

K'’T 

1 .46E+05 

8.23E+05 

1 .40E+05 

3.50E+05 

9.64E+04 

1 . lOE+01 

(continued  on  next  page) 
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TABLE  1.  (continued) 


G 

H 

I 

J 

K 

H(cm) 

1000 

1000 

1000 

1000 

1000 

SI 

0.5 

0.5 

0.5 

0.5 

0.5 

T1 ( xC) 

a. 5 

a. 5 

a. 5 

a. 5 

a. 5 

Sz ( 1 /cm) 

-^♦E-05 

-4E-05 

-4E-05 

-4E-05 

a.3aE-04 

T 2  <  xC/cm ) 

-0 . 00068 

-0 . 00068 

-0.00068 

-0.00068 

0.00068 

K ( cm/s ) 

0.01 

0.001 

0.01 

0.01 

0.01 

g 

0.56 

0.56 

0.7 

0.56 

0.56 

K ' ( cm/s ) 

o.ooa 

o.ooa 

o.ooa 

o.ooa 

o.ooa 

g ' 

0.56 

0.56 

0.56 

0.56 

0.56 

A(cm/s) 

o.oa 

o.oa 

o.oa 

0.1 

o.oa 

k 

-0.01 

-0.01 

-0.01 

-0.01 

-0.01 

h 

-0.016 

-0 . 004 

-0.004 

-0.004 

-0.004 

aS+ 

0.353 

0.874 

0.144 

0.303 

0.144 

aT+( xC) 

a.  084 

5.364 

1  .614 

a.  434 

1  .613 

aS- 

0.173 

0.794 

0.064 

o.aa3 

0.608 

aT-( xC) 

0.784 

3.904 

0.354 

1  .074 

3.973 

Ll<cm/s) 

3.9a 

1  .59 

1.59 

a.  08 

1  .59 

V  <  cm/s ) 

-1  .5 

-a.  46 

-a.  46 

-0.93 

-a.  46 

M(cm/s> 

0.015a 

0.0061 

0.0061 

0.0171 

0.0061 

L(km) 

5.15 

0.98 

4.47 

a. 15 

1  .84 

H' (cm) 

35a.  1 

430.8 

430.8 

458.9 

430.8 

S  Flux 

0.4175 

0 . 6630 

0.0837 

0.3735 

0.3989 

K"S 

4.30E+05 

1 .30E+05 

7 . 39E+04 

1 . 18E+05 

1 . lOE+05 

T  Flux 

a. 7518 

3.6443 

0.7435 

1  .8343 

1 .8831 

K"T 

5.67E+05 

1 .43E+05 

1 . 33E+05 

1 .57E+05 

1 .34E+05 

83 


